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Abstract 

In this paper we show an explicit relation between the local and the 
effective elastic properties for planar elastic structures, up to some purely 
geometric constants. The model considered is periodic structures with 
isotropic and homogeneous (or in a way uniformly shifted) local elasticity 
tensors in the linear theory of elasticity. Our study is restricted to the 
structures that have some bounded and connected periodicity cell with 
Lipschitz continuous boundary, that in addition have positive definite ef- 
fective tensors. The result is a slight generalization of previous work by 
Vigdergauz that can be found in [25] . 

1 Introduction 

In this paper we will be confined with effective elastic properties of periodic 
structures in the planar linear theory of elasticity. We show that the relation 
between the effective and local elastic properties for structures on Lipschitz 
domains is always elementary under the hypothesis of positive definite effective 
tensor and constant isotropic local tensor. This generalizes previous work by 
Vigdergauz. See [25]. A partial generalization was given in [T5] . 

The periodic structures we will consider are the planar structures with locally 
isotropic elasticity tensors which are constant throughout the domain. The 
elastic hypothesis will be the linear theory in which both the strain mapping 
and the relation between stress and strain are supposed to be linear. 

We will focus on a periodicity cell f2 of the structure to be studied. Apart 
from specifying an elasticity tensor on 17, we will make the following geometric 
assumptions on the domain. Let ft be a bounded and connected open set in R 2 
with Lipschitz continuous boundary. Moreover, we suppose that £1 is contained 
in some translation Y of the cell (0, h) x (0, 12) for some h,l2 > 0. Let be the 
periodic extension of the periodicity cell f2, that is f2 = Int Cl(£! + l\L x I2Z). 
We assume that f2 is connected and has Lipschitz continuous boundary, which 
additionally reduces the set of possible domains. 

We say that a vector u G iJ 1 (fi) 2 is quasiperiodic if u — € i/^ or (0) 2 for 
some £ e R 2x2 , where £^ er (Q) 2 denotes the closure in H 1 ^) 2 oF\C^.(Y) 2 . 
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The matrix £ is then, despite a slight risk of confusion, called the quasiperiod 
of u since the characteristic lengths l\ and I2 will be constant throughout. 

For cartesian tensor fields we will use the following notation. The set of 
symmetric n x n real matrices is identified with R™(™ -1 )/ 2 . The set of triclinic 
fourth order cartesian tensors is identified with R 6 , that is those a for which 
o-ijkl = a jiki = o-kUj for all indices. Symmetric tensors are identified with 
symmetric matrices. The elements of R 6 for which the matrix is positive definite 
will be denoted by E. The elasticity tensors will eventually be assumed to live 
in a subspace of the strongly elliptic ones £, that is those a in the Lebesgue 
space of O with values in E such that[^]£ 2 ~ £et£ on R 3 . Also, E is identified 
with its corresponding subspace of £ . An elasticity tensor a € £ is said to be 
square symmetric if aim = 02222 and aan — for all i, and it is said to bo 
isotropic if in addition aim = 01122 + 2am2- If a € £ is constant on i7, it is 
called homogeneous; for example a € E. 

We will use the linear operators on 2?'(S1) 2 obtained from the symmetric and 
antisymmetric parts of the gradient to model the linear strain and rotation of 
a displacement field u, respectively: V = e + u>. The stress a is assumed to 
be linearly related to the strain s by some a € £ via the Hooke law: a = ae, 
provided the product is defined; componentwise, 

Cjj = 22 a ijkl e kl- 
kl 

For convenience, when there is a displacement field u and an elasticity tensor 
a present, we will write just a and e for the stress <r(u) and the strain s(u), 
respectively. Otherwise, a may also denote some symbol in an equation that is 
not a priori coming from a displacement field. 

The effective elasticity tensor is defined as follows. With a € £, let a* be 
defined on R 3 by the equation 

where = a(u^) for some quasiperiodic minimizer of 3 e(u) ■ <r(u) dx 
with quasiperiod £. The above definition of a* can be justified by the Riesz 
representation theorem and homogenization theory, see for example |17) . 

Throughout, the inverse of an elasticity tensor a will be denoted by c = a^ 1 , 
and c will be called the compliance tensor. The difference between two vectors 
depending on the elasticity tensor will be denoted by the addition of a prefix 
A. For example Acim = c^n — c mi- 

The first assertion in the following theorem is the main result of this paper. 



2 The asymptotic equivalence relation ~ of which (a, b) is an element if there exists positive 
real constants Ki such that Kia < b < K2C1. 

3 The Dirichlet problem with linear boundary conditions. 
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Theorem 1. Let a € E be a constant and isotropic elasticity tensor on Q. 
Let a* be the effective elasticity tensor defined by equation |Ip with local tensor 
a. Suppose that a* is positive definite. Then there exists d £ R 6 7 which is 
independent of a, such that 

c* = c ml d + c 12 i2 e, 

where e £ R 6 is defined by e\2\2 — —%&ii22 — 1 with the rest of the compo- 
nents equal to zero. Relaxing to isotropic and smooth a £ £ , but requering that 
Acn22 = — 2AC1212 = constant, then 

Ac* = Aci2i2 e = Ac. 

A proof of Theorem [l] will be presented in Section [3j The following remarks 
concern the components of the tensor d and relations to common notation. 

Remark 1. The equation Ac* = Aci 2 i2 e in Theorem^means that for isotropic 
and not necessarily constant compliance tensor c, a shift of c such that Acn 2 2 = 
— 2Aci2i2 — constant, leads to precisely the same shift in the effective compliance 
tensor c* . This is known (see for example Muskhelishvili 's exposition on stress 
invariance in ]!&})■ See also Remark^ in this regard. 

Remark 2. The components of d in Theorem^satisfy the following inequalities: 

^1111, ^2222 > 0, 
^1111 + 2(in22 + ^2222 > 0, 
^2222 + 4(^2212 + 4g?1212 > 0, 

dim + 4dni2 + 4c?i2i2 > 0, 

This is a direct consequence of the positive definiteness of the effective tensor. 

Remark 3 (Matrix notation). In a common matri^representation of elasticity 
tensor fields, the first assertion of Theorem [7] can be stated as follows: There 
exists a symmetric matrix D with constant real entries such that the effective 
complianc^ tensor C* satisfies 

<C{ C* 2 C* 3 \ , l ,/D 1 D 2 D 3 \ / -1/2 

C* 2 C% C*\ =(- + -)[ D 2 D 4 D 5 \+-[ -1/2 
K C* 3 C* 5 C*J V A U/ \D 3 D 5 D 6 J V 1, 

where the constants K and G are local planar bulk modulus and the local shear 
modulus, respectively. Moreover, D does not depend on any of K and G. 

In the case of plane stress or plane strain, where dimension reduction or 
decoupling has preceded the planar model, the local moduli K and G above are 
related to the Young modulus E and Poisson number v as follows: 

K = — -, (plane stress) 

2(1 - v) 

K = — — -. (plane strain) 

2(l + i/)(l-2i/) ^ ' 





<^1212 


>o, 


d\\\\ - 


- 2rfii22 + ^2222 


>o, 


'2222 — 


4c?2212 + 4g?1212 


>o, 


!im — 


4dm2 + 4di212 


> 0. 



4 The so-called contracted or Voigt notation. Explicitly, C\ = cj ln , C| = cj 122 , C3 = 
2 c iii2> C4 = c 2222> @5 = ^ C 2212' ^6 = ^ C 1212' with the rest of the components implicitly 
defined by the triclinic symmetry. See Corollary [l] in Section [2] for a concrete example. 

5 The inverse C* of the effective elasticity tensor, that is effectively e' = C*tr' . 
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In either case, the shear modulus G and the corresponding three-dimensional 
bulk modulus k satisfy 



G 



2(1 + j/)' 



E 



k 



K 



G 
~3 



The Lame moduli are related by (/z, A) = (G,K — G). For definitions, we refer 
the reader to any standard reference on linear elasticity as for example JTJ Wj or 
more specifically to ' e 2£, Appendix BJ. 

We conclude this section with an example to motivate why the validity of the 
claimed result can be expected. Consider a smooth perforated periodic structure 
in R 2 with constant and isotropic local compliance tensor C and with effective 
compliance tensor C* . The classical Vigdergauz equations (see [23 Section 4]) 
can be written on the form 



where are suitable linear combinations of the Vigdergauz constants, which 
are independent of the local compliance tensor C, and Ei are constant matrices 
independent of both C and the geometry of the holes. A natural guess of the 
general form of this relation can therefore be 



for some matrix D independent of C, and some independent constant matrix 
E. An equivalent formulation, by the local isotropy, is C* ~ C\D' + CqE' , with 
matrices D' and E' corresponding to D and E, respectively, which up to scaling 
is the form we have used in the formulations above (Theorem [l] and Remark [3]). 

The result referred to in the previous paragraph is a generalization for con- 
stant isotropic tensors of the observation that if the local compliance tensor is 
changed from C to C' in such a way that C' 2 — C 2 = 2(C 6 — C 6 ) = constant, 
then the effective compliance tensor is shifted from C* to (C*)' in the very same 
way. This means that 



which agrees with equation (|2| if C is constant on fl since the above shift implies 
Cy — C\ = by isotropy. Cherkaev, Lurie, and Milton gave a proof of this result 
in [B] . As a part of the proof of Theorem [lj we provide an alternative proof of 
this. 

In the next section some direct consequences of the result will be laid out 
before we present a proof in Section [3j 

2 Square symmetry 

In this section we present a particular case of the main result (Theorem [lj in 
which the effective tensor that have additional structural symmetry. In this 
way we can relate our result to previous investigations. We remark that the 
corollary below generalizes the corresponding cited results. The references are 
included in the end of this section. We will here restrict our attention to the 



(CT.cs, cs) = c 1 (d 1 ,d 2 ,d 6 ) + c 2 (e 1 ,e 2 ,e 6 ), 



C* = C 1 D + C 2 E, 



(2) 
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effective tensors which have a square symmetric structure, that is the effective 
tensor a*,. fci satisfies a*m = 02222; an d 01112 = a 22i2 = 0- 

Knowledge about the structure of the effective tensor can sometimes be 
acquired a priori by identification of symmetry groups of the periodicity cell for 
which the inner product, corresponding to the Riesz representation of the elastic 
energy functional, is in a suitable sense invariant (see any standard reference 
on Fourier analysis and representations of groups of finite order, and for this 
particular application, see for example [SJ [Ml H3 ) • The standard example 
here is the square honeycomb, with a homogeneous and isotropic local elasticity 
tensor, for which the symmetry group is the dihedral group Dili4 of order 8. 

Corollary 1 (Square symmetry). Let the local elasticity tensor be isotropic. 
Then in terms of the planar bulk modulus K and shear modulus G, it may be 
represented^ by the matrix 

'K + G K-G s 
K-G K+G 
G) 

Suppose that the effective tensor has square symmetry, which means that it may 
be represented by a matrix of the form 



(3) 



where K* is the effective planar bulk modulus, and G* and G\ 5 are the effective 
shear moduli. Suppose that K* , G* , G\^ > 0, which is true if and only if the ma- 
trix is positive definite. Then there exist nonnegative real numbers A\, A2, A3, 
that do not depend on the local elastic properties K and G, such that 




I 

K* 


f 

_ K 


f Ai\ 




1 

*~ G 


1 

G* 


1 

~ G ' 






G, 


1 


1 

~ G ' 


M 3 ( 







Proof. Since the effective elasticity tensor, call it B* , is assumed to be square 
symmetric, it can be written on the form (|3]). The compliance tensor C* — 
(B*)^ 1 is therefore of the same form (the same entries are zero in general). The 
validity of the three equations in Corollary [l] then follows from the equation in 
Remark [3] (see also Remark [4] below) . □ 

The result of Corollary [T] is the Vigdergauz equations for the studied class 
of structures. The nonnegativity of the constants A\, A2, A3 follows from the 
assumption that the effective tensor is positive definite. We also refer the reader 
to O . To conclude this section two remarks on the geometric constants are 
included. 



3 See Remark [3] 
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Remark 4. The explicit relationship between the entries Di of the matrix D 
in Remark^ and the Vigdergauz constants A\, A 2 , A3 for square symmetric 
effective tensors, is by Corollary [JJ given by 

D l =D i = 1 -±^^, D 3 = D 5 = 0, D 2 = 1 -±^^, D e =A 3 . 

In the case of effective isotropy, Dq = A 2 , so then we have the additional de- 
pendence Dq = 2(£>i — D 2 ). 

Remark 5. We note that by Theorem^ the effective tensor a* is an elemen- 
tary function of the components of the local isotropic constant tensor a. For 
example, for a square symmetric effective tensor it follows from Corollary [JJ 
that V (K<G) K* = g(l + A u A!), V (K , G) G* = g(A 2 ,l + A 2 ), V (KG) G* 5 = 
g(Az,l + A3), where g(p,q) — (qK + pG)~ 2 (pG 2 , qK 2 ) . Thus in this case the 
effective parameters are monotonically increasing functions of the local param- 
eters K and G. In particular, K* is strictly increasing in K , and G* and G* i5 
are strictly increasing in G. For Q =£Y, the rest of the dependencies are strictly 
increasing as well, where A±, A 2 , A3 > for Q, ^ Y in Corollary^ follows for 
example from the Hashin-Shtrikman bounds. See J^j. 



3 Proof of Theorem [T] 

In this section we will give a proof of Theorem [T] First we state and prove two 
preliminary lemmas on the relation between average stresses and quasiperiods 
for the displacement fields involved in the definition of the effective tensor. 

For the proof we will make use of curves connecting opposite sides of the 
cell Q. Let 71 be a smooth curve connecting two points in O with l\ as the 
difference between the X\ components, and with equal x 2 components. We 
assume that 71 is separated from any possible hole in the global structure, that 
is dist(7i,9f2) > 0. Moreover, we suppose that the x\ component of the unit 
tangent r to 71 is always positive. Such a curve 71 exists by the connectedness of 
£1 and its periodic extension fi. Let 72 be an analogue curve in the x 2 direction, 
by interchanging all the indices. An illustration of a possible pair of curves 71 
and 72 for a periodicity cell fl is shown in Figure [l| a). 

The following lemma is a variant of the classical Cesaro formula. We will 
use the notation 

JO, ifi = j, 
I 1, otherwise. 

Lemma 1 (Quasiperiods). Let a € £ be an elasticity tensor. Letu be quasiperi- 
odic on Q, with quasiperiod £ € R 2x2 . Suppose that a and u are smooth in some 
neighborhood of 7, . Then 

-IVf-Z c f a t dx +c ■ [ d ±I*!*A x dx _<-./ d(c jqat a at ) \ 

— ^ Z^/\2J V 3 rstUst U r T ^ r 3 J Q x . ^rO,x q s rj J ^ x r ux q J, 

where c = a -1 and a is the starting point of the integration along 74. 
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Proof. Since u is quasiperiodic and smooth in a neighborhood of and the 
endpoints a and /3 of ji have equal components except for the zths, we have by 
the Newton-Leibniz formula, 



^ijU = Uj(/3) — Uj(a) = I / Vudx I = I / eefcc J + I / wdi 



By the linear relation between e and tr, and the nonsingularity of the elasticity 
tensor a, we have 



/ edxj = e jr dx r = y^ y / C J>st (T s 

■'7. / i -^7; r«t ^7i 



dx r . (4) 



By the smoothness of u, an integration by parts yields 

kr 



/ w dx J = / Wj> cfa; r = / — a r )wj r (a) — / 

r L ■ J 7,' 



dx k 



x y dx ^ 



kr,r^j *^ 



duJ r j 

— X r (IX 

ox k 



where we in the second step used that uj is periodic since u is quasiperiodic. By 
noting that 

doj rj _ de rk _ de jk 
dxk dxj dx r 



wc find 



diUrj x dxk - ( (^±_ d ^k\ x dxh 

. ^ CJX i" CJXr * 



E/ d(c r k s to- s t) d(c 3 ksto- s t) \ 
, t &r— 



El f d(CrkstO-st) , /" d(CjkstO- s t) , 
V J 9aT ^ ~~ / £b ^ : 



Thus 



,1 > / \ , V^/ /" d{c r kstO- a t) , f dCjkstO-st , 
jj (i.r j ' I 7 ( y X r OEXj - J — X r dx k 



This together with equation Q give the asserted equality. □ 

Remark 6. Note that the local elasticity tensor a in Lemma [7] need not be 
isotropic. Moreover, the lemma is true as it is stated for the analogue hypothesis 
in Ft™ for n > 2, if the factor 1/2 is replaced by l/n on the right hand side in 
the equation. 

Let U be a connected and open set in R™ with Lipschitz continuous bound- 
ary. The following version of the Green formula on ff(div, U) x iJ 1 (C7) 9 (v, ip) 
holds: 

/ v-Wipdx+ / div v ■ <p dx = (vv, cp) , (5) 
Ju Ju 
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where H(div, U) denotes the space of elements of L 2 (U) 2 with with divw in 
L 2 (U), and (•,•) denotes the pairing of H 1 / 2 (dU) and its continuous dual. 
See [11] for a construction of such a trace mapping. 




(a) (b) (c) 

Figure 1: (a) A periodicity cell fl with curves 71 and 72. (b) Translations 7^ 
and 7 2 of 71 and 72, respectively, (c) The domain with outer boundary T 
corresponding to the curves in (a) and (b); fi' has no cusps. 




Figure 2: An interior part of a periodically perforated structure. 

We will make use of the following curves in the extended domain. Let "/[ be 
the translate of 71 such that the left endpoint of r y' 1 is the common point of ji 
and 72. Then let 7 2 be the translate of 72 such that its lower endpoint is the 
left endpoint of j[ . Let Y' be a translate of Y such that and j' 2 are the lower 
and left sides of the boundary of Y' n A, respectively. Let O' be the intersection 
of fi and the region of R 2 bounded by 

T := 7 J U 7 2 U(0,Z 2 )+ 7 iU (?!,()) +7 2 - 

Illustrations of such curves and domains are shown in Figure [3a)-(c) for the 
periodically perforated structure shown in Figure [2j 



8 



Lemma 2 (Average stress). Let the local elasticity tensor in £ be smooth on f2 
and Y -periodic. Suppose that u is a quasiperiodic minimizer of J Q e{u)-a{u) dx. 
Then for i ^ j, 

/ an dx — li I [av)i dx, / <Xy dx = lj / (av)i dx. 

Jn Jfj Jn Jfi 

Proof. Since u is quasiperiodic and minimizes the elastic energy, u is the unique 
solution, up to translation, to the equation 

diver = in V'(h) 2 , 
ais = in H- 1/2 (dn) 2 . 

Since a and diver have components in L 2 (fl'), it follows from the Green 
formula ([5| that 

(0-11,012) • Vipdx+ / div(crii,cri 2 ) fdx = ((av)i,(f), 
Jn' 

for any ip £ _ff 1 (f2'). By the regularity of u and the vanishing of av on <9f2, we 
have that av G L 2 (dfl') 2 . See [5]. Hence 

{(crv)\,<p) = / [crv)\^pdx = / (au)i(pdx 
Jdci' Jr 

(crv)iip dx + / (av)i(pdx+ / {o~v)i <p dx + / (av)\ipdx 

-/(0,i 2 )+7l ^(d.0)+72 

(cr!/)i (y>(as) - (p(xi,x 2 - h)) dx+ (av) l (ip(x) - (p(x 1 - h,x 2 )) dx, 

7i "'72 

where we in the last step made a change of variables and used the periodicity 
of a and the opposite signs of the outward unit normals. 

Let if = x\. Then tp(xi, x 2 — 12) — fip) an d vC 3 -) — fi x i ~ h, x 2 ) = h- Thus 



((au)\, ip) = h / {av)\dx — l\ I (av)idx, 

where we in the last step used the periodicity of a. Since Vip = (1,0) and 
div cr = 0, we have by periodicity and the Green formula that 

/ o"ii dx = / en dx = l\ I {<7v)idx. (6) 
Jn Jn' Jj 2 

With tp = x 2 , we obtain in the same way, 

/ ai 2 dx = l 2 / (av)idx. (7) 
Jn J ^ 

By the symmetry with respect to the indices, the equations corresponding 
to the equations ^ and §f§ are 

/ a 22 dx = l 2 / {av) 2 dx, / a 12 dx = h / i?v) 2 dx, 

J n ^71 J n J 72 

which completes the proof. □ 
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Now we shall specialize Lemma [T] to the case of constant and isotropic local 
elasticity tensors. We will us^] the following notation for scalar (j> and vector- 
valued ip: 

_ ( d<f> d(j) \ _ d(fi 2 dipi 



curlt/) = —,-— , curl<^ = 

\OX2 OX\ ) OX\ OX2 

which is the reason for the slight antisymmetry in the following formulas, with 
respect to the indices. 

Lemma 3 (Quasiperiods for constant isotropic tensors). Under the hypotheses 
of Lemma [7J suppose that a is constant and isotropic, £ is symmetric, and that 
div a — on 71 U 72 • Then 

£11 = ~ ^nsT I CT 22 dx + / ((Tr o, 0) - x 2 curl Tr a) dx, 

\ Y \ Jn n J 7l 

£22 — - 2 ,' 2 ' 2 / <7xi dx + 01111 / ((0, Tr a) + x x curl Trc) dx, 
l y l '2 J 7, 



_ 2C1212 

62 -~jyr 



2Z 2 

Proof. By Lemma [I] we have 



/ CT12 dx H — / ((Tr cr, 0) + xi curl Tr er) dx 
/ ((0, Tr a) — X2 curl Trcr) dx. 



Cllll 



CiiA = ^ / ci rst o- s tdx r + ( f c 2qs S^-X2dx q - I 

rst J li gat ^•'71 A 



5(7 



Clqst : Q^-X2 dx q 



Since a is isotropic, 

Cllll = C2222, Clin = C1122 + 2C1212, C1112 = C2212 = 0. 

By Lemma [2] we have 

t - f o~22 dx == f (av) 2 dx^ J er 2 2 dxi - j o- 12 dx 2 , 



— / o\\dx= I (av)idx= I un dx2 — I a^dxi. 

'\ J £1 J 72 </ 72 72 



7 This double notation (u> and curl) for the so-called rotation of vector fields with values in 
R, 2 is motivated by the standard notation, since we will here not be concerned with specific 
components. 
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The representation of £n is obtained by the following calculation: 

£llZl = Cuil/ O-n dx X + C U2 2 / CT 22 dx 1 + 2C1212 / CT 12 efo 2 
•/71 "'71 •'71 

+ 2c m2 / 17^2 ^1 + C1122 / ^ ll 2 ; 2 ^2 + cnii / ^r^2 ^2 
J lt oxi J lt dxi J lt dxi 

f 9cr H J A 9cr 22 j n [ 5(J 12 j 

-cuu / — — x 2 dxi-cn 22 / — — x 2 dx 1 -2c 12 i2 / ^ — x 2 dx 2 
7 7l cte2 7 7l «r 2 J 7l dx 2 

= cuu / cru da:i + (cim - 201212) / cr 22 cfai + 2c m2 / a 12 dx 2 

+ 2c m2 / ^p^-x-z dxi + (cim - 2C1212) / ^-^£2 ^2 
J 71 O^i J 71 OTl 

/" (9cr 22 /" (9(7ii 

+ C1111 / — x 2 dx 2 - cim / — — x 2 c(a;i 
7 71 oxi 7 71 cte 2 

- (cim - 2ci2ia) / ^r^a;2 dxi - 2c m2 / ^p^-x 2 dx 2 
J 71 ox 2 J 71 dx 2 

f rp j . /" 9Tra /■ dTra 

= c U n / Trcrdxi +C1111 / — x 2 dx 2 -ci in / — x 2 efei 

7 71 J 7l OXx J yi dx 2 

-201212(7- / a 22dx~- ( (divcr) 2 a; 2 o?.xi + / (div a)ix 2 dx 2 ) 

\'2 Jsi V 7l 7 7l / 

= cim / ((Tr <T, 0) — x 2 curl Tr ct) efe / a 22 dx, 

Jit '2 7n 

where we in the last step used div a — 0. By the same procedure, we find 

£22^2 = cim / ((0, Tr a) + xi curl Tr a) dx - 2ci212 f a n dx, 
J~i2 n Jo. 

from which the representation of £ 22 follows. 

We turn to £i 2 . Let a denote the common point of 71 and 72. By Lemma [2] 
and the periodicity we have 



— / 012 dx = / (av) 1 dx— I o\ 2 dx\— f o\\dx 2 , 
2 Jn J~f{ J-/[ Jj[ 

- / o 12 dx= I (av) 2 dx= / o\ 2 dx 2 - / a 22 dx\. 

'1 ill J 7 2 «/ 7 2 •'72 
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By periodicity it follows from Lemma [T] that 
(£12 +W2i(a))Zi 

= 2ci2i2 / 012 dxx + (cim - 2ci2ia) / 0"n dx 2 + cim / a 22 dx 2 
+ C1111 / ^p^xi dxi + (cim - 2ci2i 2 ) / ^^xidxi 

J-yJ C/X2 J-yJ CX2 

f 9ai2 f doi2 

+ 2C1212 / -5 Xi dx 2 - 2C1212 / -5 — SCl ctai 

9x 2 J 7 j dxi 

- (cim - 2ci2i 2 ) J ~^r x i dx2 - cim ^ d:E 2 

9 Tr a , f d Tr a 



= cim / Trcrefe 2 + cxxxi / — x x dxi - cnu / — xi cix 2 

7 7 ; J 7 ; ox 2 Jy ± 0x1 

+ 2C1212 J a 12 dx - J (div a) 2 x 1 dx 1 + J (div (t)i x : dx 2 ^ 
= Cnxi / ((0, Trcr) + Xi curlTr <r) dx H ^ 212 / a 12 dx. 



Similarly, 

(61 + ui 2 (a))l 2 — cnn / ((Tr cr, 0) — X2 curl Tr a) dx H f a 12 dx. 

Jy 2 n Jn 

Since £ € R 3 and w is antisymmetric, we have 

2£i2 = I ((Tr cr, 0) + xx curl Tr cr) dx 

k J 7l 

+ / ((0, Trcr) — x 2 curl Trcr) dx H r^^- / a 12 dx, 

'2 J 7 ^ |J| Jo 

which completes the proof. □ 

For the uniform shifts we will need the following version of the above lemma. 

Lemma 4 (Quasiperiods for uniformly shifted isotropic tensors). Under the 
hypotheses of Lemma [7J suppose that a is isotropic, £ is symmetric, and that 
diver = on 7! U 72. Assume that c is shifted in such a way that Acn 2 2 = 
— 2Aci2i2 = constant. Then 

2Ac 1212 f 
Atxx = -^ vr J n ^dx, 

2Ac 1212 f 
2Ac 1212 f 

Proof. By considering the differences A^j and carrying out the same calcula- 
tions as in the proof of Lemma [3j the result is obtained by noting that the 
supposed shift for isotropic c means, in particular, that Acim =0. □ 
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For presentation purpose, we include a fundamental lemma which is some- 
times attributed^ to Michell and it is concerned with the dependence of the 
stress on the elasticity tensor for the Neumann problem. See [15]. This is the 
reason why we restrict our analysis to two dimensions. The version we are 
interested in can be stated as follows. 

Lemma 5 (Michell lemma). Let the elasticity tensor a £ £ be smooth and 
isotropic on a bounded and connected Lipschitz domain fl in R 2 . Then the 
stress a does not depend on the elasticity tensor in the pure Neumann problem 
on with (ov, l)r = for all connected components T of dVL, 

(i) if a is constant, or 
(ii) i/Acn22 = — 2Aci2i2 = constant. 

In (i) the condition on the data is necessary if A(a/ami) ^ 0. On the other 
hand, if a is constant and A(a/aim) = 0, the condition on the data can be 
relaxed to just equilibrium, (av, I) en = 0, and the stress does not depend on the 
specific components of the elasticity tensor. 

Remark 7. The property A(a/ann) 7^ means that two constant elasticity 
tensors are not scalar multiples of each other, which is the case if and only if 
A(cii22/cini) 7^ 0. See also Remark^ 

Before giving a proof of the lemma we make some additional technical re- 
marks. The boundary data in Lemma [5] should be understood to belong to 
i/~V 2 (T\) anc j j g na turally assumed to satisfy J ^2 i (o'i / ,i)r i = in either case. 
Here (•, -)r 4 denotes the pairing of H 1 / 2 ^^ 2 and its continuous dual, where Fj 
are the connected components of <9f2 with Tq being the outer part. 

By the supposed connectedness of ft, an element v g L 2 (£!) 2 satisfies div v = 
and (vu, l) H i/2^ i ) = if and only if there exists a unique <fi £ H 1 (fl)/'R such 
that v — curl (j). Moreover, the following Hclmholtz type decomposition holds: 

L 2 {Q) 2 = curl{0 e H\n) : 0, ro = 0, <j>\ Ti = constant} © V-ff 1 ^) (8) 

where the sum is direct. Proofs of these characterizations can be found in (TTJ 
Chapter 2]. 

We will also make use of a necessary and sufficient compatibility condition 
on symmetric Vij £ L 2 (fl) that guarantees the existence of some displacement 
field u £ i? 1 (r2) 2 such that v — e(u). We will call the condition the Donatj^] 
lemma and it can be stated as 



for all symmetric s with components, by density, in C5 (f2) and such that div s = 



Remark 8. Note that f2 is not required to be simply connected in the Donati 
lemma, only connected, bounded, and open in R 2 , and with Lipschitz continuous 




(9) 



0. See 




See 1161 Chapter 5] for some historical remarks. 
Sometimes attributed to Donati in this context. 
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Proof of Lemma^ By an argument by Frechet and Riesz (PHHH]), there exists 
a unique u £ H 1 ^) 2 / kei(e) such that 

diva = in V'(fl) 2 , 
uv = gi iaH- 1 / 2 ^) 2 , 

for any such gi satisfying (g il l)Y i = 0. Let a be the solution to the above 
problem for some given gi with constant and isotropic elasticity tensor on Q. 
Then a has components in L 2 (Vt) and thus there exists <f> £ i? 1 (f2) 2 such that 



dx-2 



0~21 



dxi 



0~12 



d(j)i _ d(j> 2 
dxi dx 2 ' 



(10) 



where the last equation comes from the symmetry of a. The condition (J9j) can 
be written as follows in terms of 4>, for any symmetric s with components in 
C§°(£1) satisfying divs = 0, 



Indeed, 



= — / SijSij dx — cim / curl <fi Tr s dx. 



^ CijkiOkiSij 

ijkl 

(cilllCTll + Cii22Cr22)sil + (C1122C11 + C1H1CT22)S22 + 4ci 2 12Cl2Sl2 



Cllll 



+ 2C12 



uepi _ uq>2 
8x2 dxi 



dxi 



Sll 



(s n + s 22 ) 



oq>2 oepi 
~ — S12 - 3 — S12 
0x2 0x1 



dx' 



-S22 



where we in the last step used the equations in (10) and the isotropy cim = 
C1122 + 2ci2i2- Moreover, 



and 



dxi 



dxi 



Sll 



S12 



dx 2 



dx 2 



S12 dx 



dsn dsi 



dxi dx 2 



S22 dx 



dsi2 
dxi 



9S22\ 

dx 2 ) 



dx = 0, 



dx = 0, 



where both integrals vanish because div s = 0. Thus 



? i £ijSij dx = 01111 J n {w 2 - (sn + S22) dx ' 



as claimed. Hence J n curl <p Tr s dx — 0, which by the sufficiency of equation (|9j) 
in the Donati lemma, guarantees the existence of some v e H 1 ^) 2 such that 
the given a comes from v. Since this condition does not depend on the elas- 
ticity tensor, the displacement field v for any supposed elasticity tensor can be 
recovered by using the same <f>, which proves the first assertion (i). 
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We turn to (ii) and note that by isotropy Acim = AC1122 + 2AC1212 = 0. 
By the calculation in (i), we thus have 



( ^02 9^2 94>i 84>i 

ij = 2AC 12 12 q S U + - S12 - ~ S12 - ^ «22 

V oxi 0x2 ox\ axi 

%j 



Since AC1212 is constant, we conclude that 

■■ij^ij d>X 0, 



A E 



which shows the claimed independence. 

We consider now the case when the condition on the data is relaxed, but 
a is constant and A (01122/01111) = 0. Since J2i(crv, l)r, = 0, the stress field a 
exists with components in L 2 (f2) and by the Helmholtz decomposition Q there 
exist <fi, q G i/ 1 (r2) 2 such that 

_ _ \ / 94>i d<t>i \ 

Til 0-12 \ _ I 0x 2 dxi I , y„ 

By the necessity of equation we have for any s as above, that 

/ £ijSij dx = Cuii / curl Tr s dx — ci 122 / divgTrsdx. (11) 
' in Jn Jn 



E 



Provided (erf, l)r ; 7^ for some i, the equation (11) holds independently of 
the elasticity tensor if and only if A(cu22/cuu) = 0, since dive/ = gives 
self-equilibriated holes, (eri/, l) r . = 0. Thus, in general, then and only then 
there exists a displacement field v € iif 1 (SI) 2 coming from a, by the Donati 
lemma (|9|. □ 

Remark 9. The property (ii) in Lemma^means invariance of the stress in the 
studied boundary value problem under the family of maps 

1 1 \ ( 1 1 



where K and G denote the planar bulk modulus and the shear modulus. This 
particular invariance was investigated in J^. We refer the reader to 1181 \iyf for 
a thorough study of invariance of solutions in planar elasticity. Moreover, the 
fraction C1122/C1111 is an elementary function of K/G, as well as of the Lame 
ratio A//i, which are canonical parameters in planar linear elasticity. See 11 61 

W)- 

Now we are in a position to prove the main result. 

Proof of Theorem^R By the hypothesis of positive definiteness on R 3 of a* de- 
fined by equation (fTl) , a* is nonsingular since by the definition it is symmetric. 
Hence, for every average stress there exist a unique stress field a that, by the 
Michell lemma (Lemma [5]) , does not depend on the isotropic local elasticity 
tensor a £ E. For the boundary of any hole in the global problem is stress-free 
and hence the average normal stress on the outer part Tq of the boundary of f2 
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vanishes by the Green formula. The corresponding quasiperiodic displacement 
field u is unique up to translation. Thus the quasiperiod £ G R 3 of u is uniquely 
determined by a and a. 

Let (a) denote the average stress |y | — 1 J Q (an, er 22 , 012) dx, which we will 
let vary over the canonical basis /j of R 3 . With (a) = fi, equation gives 
c* m = ((a*) _1 )ini = £n- Since a is isotropic, £ G R 3 , and diva vanishes, 
there exists by Lemma [3] a real constant dim, which depends only on a and 
hence not on a, such that c ml = cim dim. Moreover, by the same equation, 
there exists constants dn 22 and dm2 such that c* 122 = — 2ci2i2+cmi dim and 
c iii2 = c nn ^1112- With (cr) = / 2 , we similarly find constants e? 2222 and <i 22 i 2 
such that c 2222 = cnn rf 2 222 and c 2212 = c m id 22 i 2 . Finally, with (cr) = / 3 , 
we find a constant dm2 such that c m2 — C1212 + cim g?i 2 i 2 . Let d G R 6 with 
components dijki defined as above and extended by zero, up to the triclinic 
symmetry. Let e G R 6 be defined by ei 2 i2 = — 2en22 = 1, with the other 
components set to zero. Then c* = cim d + cyin e> where d and e do not 
depend on a by construction. 

We argue in the same way for the case when a is changed in such a way that 
Acii 2 2 = — 2AC1212 = constant. First, we have Acim = by isotropy, and 
by the assumption of uniform shift of the not necessarily constant compliance 
tensor c, the stress is invariant by (ii) in the Michell lemma. By Lemma |4j we 
have Ac mi = A£n = 0, Ac^ 122 = -2Ac 12 i 2 = Acn 2 2, and Ac m2 = 0, for 
(a) = f\. For (cr) = f2, wc find Ac 2222 = and Ac 2212 = 0. Finally, we find 
^ c i2i2 — A c i2i2 for (a) = fa. Hence, Ac* = Aci 2 i 2 e = Ac. □ 

4 Representations of geometric moduli 

For computational purpose it may be worthwhile to document some representa- 
tions of the geometric constants diju . In the proof of Theorem [l] the constants 
were represented in terms of average values of the stress field for some given 
global average stress. The reason for this was that we were interested in find- 
ing constants that were indepentent of the local elastic moduli by virtue of the 
Michell lemma. If we want to calculate the geometric constants however, it may 
be more convenient to start by specifying the quasiperiod and then calculate 
the average stress. This is because a numerical method, as for example a finite 
element method, would then be completely explicit. For this purpose we prove 
the following lemma, in which the geometric tensor dijki is represented by a 
symmetric matrix D in R 6 as described in Remark [3J Here we will use the 
notation (/) = t^t J n f dx to denote the weighted average value of a function 
over the periodicity cell fl. 

Lemma 6. Let the local elasticity tensor be the diagonal element in E with 
K = G = 1/2. Let o~ a , er b , o~ c be the stress fields corresponding to the respective 
£ Q i £, b > £ c quasiperiodic minimizers of J^e-adx, where £ a , £ & , £ c is set to be the 
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canonical basis of R 3 . Then 

D 1 =f,- 1 ({a\ 2 ) 2 -{at 2 ){4 2 )), 

D 2 = ^{{a a 22 ){al 2 ) (a a 12 ){a b 12 ) (a a n ) (<J C 12 ) (a b 22 ) 

+ Ki)(^ 2 > 2 + K 2 ) 2 K 2 > + (^ 2 >K 2 > 2 - 2K 2 )K 2 )(^ 2 )), 
A3 = /i- 1 (K 2 )(^ 2 )-(^ 2 )K 2 )), 
^ = /i- 1 (K 2 ) 2 -K 2 )K 1 )), 

^5 = M- 1 (^i2>Ki)-K 2 }K 2 » ! 

As = m- 1 ^} 2 - - 2K 1 )K 2 )( ( r 2 b 2 ) 

- 2K 1 )(a? 2 ) 2 - 2(^ 2 ) 2 K 2 ) - 2(a b 2 )K 2 ) 2 +4(a? 2 )K 2 )(^ 2 )), 
where fi is the nonzero real number defined by the equation 

Proof. For £ £ R 3 , the effective relation Q can by Theorem [l] and Remark [3J 
be written as 

(011, 0-22,0-12) dx = A* (£u, £22,2^12), 
where A* is the inverse of the nonsingular matrix 

/l 1\M ft M i/0 -1/2 0\ 

Mis a aJ +5 U /2 :] 

By letting £ vary over the canonical basis of R 3 we obtain a system of linear 
equations in the components of D: 

B(D 1 , . . . ,D 6 ) = J^o-^a^o-^, a b 2 ,a c 12 ) dx 

This system is uniquely solvable with solution D given explicitly by the equa- 
tions in the statement of the lemma. □ 

Remark 10. The reason for the supposition that K = G = 1/2 in Lemma^ 
is only notational convenience. The same argument gives the corresponding 
representation for arbitrary K, G > 0. Also, for numerical calculations it may 
be convenient to have a local elasticity tensor which is diagonal, which here 
corresponds to zero Poisson ratio in a planar model. 

We document the representations found in the proof of Theorem [T] in the 
form of a lemma. 

Lemma 7. Let K and G be the local elastic moduli. Let a a , a b , a c and £ a , £ b , 
£ c be the stress fields and quasiperiods corresponding to the quasiperiodic mini- 
mizers of J n e ■ a dx with average stress (a) equal to (1,0,0), (0,1,0), (0,0,1), 
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respectively. Then 

D l = ^-{ ((Trcr a ,0) -x 2 curlTrcr a )da;, 

\- I ((0,Tra a ) -xiCurlTrcr a )dx = -,- / ((Trer b , 0) - x 2 curlTra b )da:, 

^/ ((Trer a ,0) +x 1 cmiTro- a )(fa+ ~ [ ((0, Tr cr a ) - x 2 curl Tr <j a )dx 
h J 7l 4/ 2 J 72 

\- { ((Trcr c ,0) -x 2 curl Trcr c ) dx, 

4/1 y 71 

— / ((0,Tra b ) + xiCurlTrcr b )cfe, 

^- / ((0,Tra c ) + x x curlTr cr c ) cZz, 

D S = Tr I ((Tra c ,0) +xiCurlTrcr c )dx+ -,- / ((0, Tr cr c ) - x 2 curl Trcr c ) dx. 
4ii J 7 , 4i 2 J„ 



We conclude this section with a remark on a representation that can be used 
for numerical calculations. We will write the geometric constants in terms of 
the average elastic energies of the solutions to the cell problems. Consider the 
following cell problems with denoting the "average strain" : 

e = (1,0,0), e = (0,1,0), e = (0,0,1), 
e = (1,1,0), e = (0,1,1), e 6 = (1,0,1). 

Let u l denote the solution to equation ^ with £ — and let B* denote the 
effective elasticity tensor in the notation of Remark [3j First, with 



D 3 = 
3 Hi 



D * = 4l 



D 5 = 
5 4Z 




e{u l ) ■ Be(u j )dx, 



£11 


£12 


£13 


£12 


£22 


£23 


£l3 


£23 


£33 



we have 



By the equation in Remark[3j the entries of the matrix D can then be calculated. 
Here, only the the diagonal entries of the matrix B* were expressed as energies 
of the solutions to the cell problems. We observe that, by the definition of the 
effective elasticity tensor, 

on — -fc»i, £22 — i>4, £-33 — ±> 6 , 

044 — n 1 + La 2 + -D4 , 055 — i) 4 + zi> 5 + ±> 6 , Cqq — n 1 + Zn 3 + is 6 . 

By linearity we can find u 4 , u 5 , and w 6 , from u 1 , u 2 , and it 3 . Moreover, it follows 
that 

D * c D * £44 — £11 ~ £22 D * ^66 — £11 — £33 

B 1 =tu, B 2 - , B 3 



2 ' J 2 

^55 — ^22 — <^33 

2 



i>4 — ^22, -D5 — ^ , £% — ^33, 
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which is now explicitly written in terms of the average elastic energies En /2 of 
the solutions u % to the cell problems. Finally, the representation obtained for 
the geometric elastic modulus D is therefore 



D 



KG 



K + G 



2ei e 4 - ei - £2 

£4 - £1 - £2 2£ 2 
\ e 6 — e l — £3 e 5 ~ e 2 — e 3 




where ej = £i»/2 denotes the average elastic energy for the solution u l with 
quasiperiod £*. 



5 Numerical example 

Consider a periodic symmetric material structure consisting of one isotropic 
material of Young modulus E = 1 with Poisson ratio v (for example a plane 
strain model), which is perforated in a periodic manner. In Figure |3] one period 
(cell) of this material structure is illustrated. 

Due to the mirror symmetry in the cell and the local elasticity tensor, the 
effective elasticity tensor orthotropic, that is the effective stiffness matrix 
B* is of the form 



B* 



We may, again by symmetrjj^] solve the cell-problem on one half of the cell (see 
Figure [SI, and we set Y = (0, l) 2 . 

To find the parameter B\ we use the following Dirichlet boundary conditions: 

£ = ui(0,x 2 ) = 1*2(2:1,0) = u 2 (xi, 1) = 0, u x (l, x 2 ) = 1, 

which, as indicated, correspond to £ = i^ 1 = (1, 0, 0) as in the end of the previous 
section. That means that B\ is equal to two times the average elastic energy. 
Similarly, to find the other parameters we use the following boundary conditions: 




u 2 (xx,0) = ui(0,x 2 ) = ui(l,x 2 ) = 0, u 2 (xi, 1) = 1, 

m 2 (0, x 2 ) = ui(xi,0) = 0, u 2 (l, x 2 ) = ui(xi, 1) = 1/2, 

ui(0,x 2 ) = "2(2:1,0) = 0, ui(l,x 2 ) = u 2 (xi,l) = 1. 



After the energies £1, . . . , £4 have been calculated, we can use the formula for D 
given at the very end of the previous section. 

The numerical solutions are obtained by using the finite element software 
ANSYSp^jwith the plane strain model. This means that (E,u) and (K,G) are 
here related as described in Remark [31 We have used a 2D 8-node structural 



10 See Section pi and the references given there. 

llr The reason for not choosing a smaller part of the domain is just for notational convenience. 
Actually, all computations can be carried out on the restriction to one-fourth of the cell shown 
in Figure [3] for example (0, 1) X (0, 1/2). This is because the order of the symmetry group 
Dih2 is 4. Moreover, the group is commutative so the equations will also then decouple. 

12 Version 12.2. 
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solid element called planel83 and the element mesh consists of 9936 elements 
and 30235 nodes, as seen in Figure [4] By finding the energy e, in the four 
cases described above we obtain the numerical values for the case E = 1 (for 
the general case one can just multiply the values of B* given below by E), 
and these are presented in Table [T] From the untruncated data corresponding 
to Table [T] we find the numerical values of the parameters D, by using the last 
equation in the previous section. These values are presented in Table [2] We 
observe that the numerical values of these parameters are almost independent 
of the Poisson ratio, in agreement with our theoretical results. 



V 


10" 6 


10" 4 


0.1 


0.2 


0.3 


0.4 


0.49 


B\ 


0.756 


0.756 


0.776 


0.835 


0.970 


1.324 


2.716 


B% 


0.806 


0.806 


0.827 


0.890 


1.034 


1.412 


2.894 


B* 2 


0.037 


0.037 


0.107 


0.209 


0.382 


0.775 


2.233 


B* 6 


0.309 


0.309 


0.292 


0.279 


0.268 


0.260 


0.255 



Table 1: The numerical values (truncated) of the entries of B* obtained from 
ANSYS. 



V 


10" 6 


10" 4 


0.1 


0.2 


0.3 


0.4 


0.49 




0.33123 


0.33125 


0.33125 


0.33123 


0.33123 


0.33125 


0.33120 


D 2 


0.23466 


0.23466 


0.23466 


0.23466 


0.23466 


0.23464 


0.23468 


D 4 


0.31078 


0.31078 


0.31078 


0.31078 


0.31078 


0.31080 


0.31078 


D 6 


0.30835 


0.30835 


0.30835 


0.30835 


0.30833 


0.30833 


0.30835 



Table 2: The numerical values (rounded) of the entries of D for the various 
Poisson ratios obtained from the untruncated data corresponding to Table [T] 
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Figure 3: The cell of periodicity is a 2 x 1 rectangle with a hole of radius 1/4 
in the center of the cell. 
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Figure 4: One half of the cell. 
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